Fluorescence quenching near small metal nanoparticles 
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We develop a microscopic model for fluorescence of a molecule (or semiconductor quantum dot) 
near a small metal nanoparticle. When a molecule is situated close to metal surface, its fluorescence 
is quenched due to energy transfer to the metal. We perform quantum-mechanical calculations of 
energy transfer rates for nanometer-sized An nanoparticles and flnd that non-local and quantum-size 
effects significantly enhance dissipation in metal as compared to those predicted by semiclassical 
electromagnetic models. However, the dependence of transfer rates on molecule's distance to metal 
nanoparticle surface, d, is significantly weaker than the d~^ behavior for fiat metal surface with a 
sharp boundary predicted by previous calculations within random phase approximation. 



in 



CD 



a 



o 

> 

in 

\o 
^^ 

rn 

O 



X 



I. INTRODUCTION 



During past decade, fluorescence of molecules placed 
near a metal nanostructure supporting surface plasmon 
(SP) has attracted renewed interest due to biosensing 
applicationSi^"'^ While early studies mainly focused on 
fluorescence near rough metal films;^^— more recent ad- 
vances in near-field optics and in chemical control of 
molecule-nanostructure complexes spurred a number of 
fluorescence measurements from dye molecules and semi- 
conductor quantum dots (QDs) linked to metal nanopar- 
ticles (NP)ji2r— Fluorescence quantum efficiency of such 
conjugate systems is determined by competition between 
strong local field enhancement of molecule's dipole mo- 
ment near SP resonance^ and, hence, of radiative de- 
cay rate F^, and nonradiative energy transfer from ex- 
cited dipole to optically-inactive electronic excitations 
in the metal characterized by nonradiative decay rate 
Tnr'^'^ The interplay between these two mechanisms 
depends on separation, d, of the emitter from the metal 
surface i^i^ As a result, fluorescence is enhanced at some 
optimal d and is quenched close to NP due to the sup- 
pression of quantum efficiency, Q = Tr/ (F^ +F„r), by 
non-radiative processes. Both enhancement and quench- 
ing were widely observed in fluorescence experiments on 
Au and Ag nanoparticles.—"— 

When a dipole with excitation energy uj is situated 
near a spherical NP of radius R its energy is transferred 
to via Coulomb interactions to electronic excitations in 
NP (single-particle or collective) with angular momen- 
tum I. Within classical approach, ^^''''^ these are collec- 
tive modes described by resonances of NP Z-pole dynamic 
polarizabilities ai{uj) expressed in terms of local dielec- 
tric function e{uj,r), i.e., bulk metal dielectric function 
e{uj) — e'{uj) + ie"{ijj) ioi r < R and dielectric constant 
of outside medium £„ for r > R. The energy dissipation 
in NP is then determined by e"{uj) which comes mainly 
from electron-phonon interactions and, for high w, from 
optical transitions between d-band and sp-band. The 
distance dependence of energy transfer rate F„r stems 
from electromagnetic interaction of molecular dipole with 
its image that, at small d (i.e., d/R <^ 1), yields^i^ 



Tnr oc d~^ , similar to flat metal surface;^ Recent flu- 
orescence measurements from molecules near (but not 
very close to) relatively large (tens of nm) NPs^^"— in- 
dicated a good agreement with classical electromagnetic 
mo dels <^""— 

With decreasing distance between dipole and metal 
surface, the energy transfer rate gets significantly mod- 
ified due to spatial dispersion of e. For d < vp/^, 
where vp is Fermi velocity, the high Fourier harmon- 
ics of dipole's Coulomb potential are no longer screened 
by conduction electrons in the metal so that the phase 
space restriction on excitation of electron-hole pairs over 
Fermi level is liftedi^ For small d, electron-hole gener- 
ation becomes the dominant energy transfer mechanism 
that leads to a sharper distance dependence of energy 
transfer rate, F„r oc d"^; the precise coefficient calcu- 
lated within random phase approximation (RPA) differs 
in the literature42iii 

While fluorescence quenching from a molecule near a 
relatively large NP is essentially similar to that near flat 
metal surface,'^®''*- the energy transfer to small NPs is 
significantly modified due to quantum-size effects. Flu- 
orescence near nanometer-sized Au NPs with diameter 
as small as 1.4 nm has attracted significant interest 
due to their large surface to volume ratio and, hence, 
lower Ohmic losses, which is important for biosensing 
applications i^i^ For such NPs, a spillover of a sub- 
stantial fraction of electron density beyond NP clas- 
sical radius can significantly affect spatial dispersion 
of metal dielectric function"*^""*^ and hence modify the 
dipole Coulomb potential near NP surface. The presence 
of strongly localized d-band with abrupt nearly classical 
density profile gives rise to a surface layer where Coulomb 
potential is not screened by interband transitions.—"— 
In small NPs, the size-dependent Landau damping of 
plasmons excited by the molecular dipole can increase 
the rate of energy dissipation.^^ Incorporation of all the 
above effects in a single model requires a microscopic ap- 
proach of energy transfer rates near small NPs. 

In this paper, we present such a microscopic model, 
based on time-dependent local density approximation 
(TDLDA), for calculation of radiative and nonradiative 
decay rates of a molecular dipole near a small noble metal 



NP. Here, the molecule is represented by a point-like 
dipole while we focus on non-local and quantum-size ef- 
fects in dielectric function of nanometer-sized metal NP. 
To this end, we adopt quantum three-region model we de- 
veloped recentlji^i^ that accounts for different density 
profiles of sp-band and d-band electrona^"— in small no- 
ble metal NP embedded in a dielectric medium. We find 
that calculated radiative decay rates, F^, are not signifi- 
cantly affected by non-local effects and are comparable to 
those calculated using semiclassical approache a -^''^^ with 
quantum corrections. At the same time, the energy trans- 
fer rate Tnr is significantly enhanced due to the gener- 
ation of electron-hole pairs in metal at small molecule- 
NP separations. We find, however, that quantum-size 
and non-local effects strongly affect the distance de- 
pendence of Tnr close to NP surface and that d~^ be- 
havior no longer holds; namely, the energy transfer 
rate is significantly slower than that predicted by RPA 
calculations 4Mi 

The paper is organized as follows. In Section [ll] we 
set up the expressions for decay rates to be evaluated 
and derive the classical limit. In Section IIIII we outline 
the quantum three-region model used for evaluation of 
NP polarizabilities. In Section IIVI our numerical results 
are presented and discussed, and SectionlVlconcludes the 
paper. 



II. RADIATION OF A DIPOLE NEAR 
SPHERICAL METAL NANOPARTICLE 

Consider a molecule with a dipole moment fi located 
at position r measured from the center of a metal NP of 
radius R placed in a medium with dielectric constant e™ . 
Its decay rate is given by a general expression^ 



second term in Eq. ([3]) contains NP Z-pole dynamic po- 
larizability, a;(w), given by 
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where Sn^'"' (r) is the induced density due to potential r' 
(hereafter we suppress w-dependence to alleviate expres- 
sions). From Eqs. ([T]) and ([3|), F„r can be calculated for 
arbitrary dipole orientation with respect to the NP sur- 
face; for perpendicular and parallel orientations, it takes 
the form^sdo 
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where r° = l^^fc^ is radiative decay rate of an iso- 
lated dipole. The radiative decay rates in the presence of 
NP are similarly obtained using Green dyadic's far field 
asymptotics as^^iiS, 
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and they depend on NP dipole polarizability, ai(w), that 
peaks at surface plasmon frequency, u/gp- Note that the 
expression ((4]) for NP polarizability that determines de- 
cay rates ([5]) and ^ is valid as long as there is no di- 
rect overlap between molecule orbitals and electron states 
in the metal j^ In classical approximation of a sharp 
NP boundary at R, the induced density Sni{r) peaks at 
r — R, yielding standard expression for NP polarizability. 
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r = 47rfc2lm/x-G(r,r) • /x. 
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where G'pi,(r,r') is the electric field Green dyadic in the 
presence of NP and k is the wave vector of light. The full 
decay rate Eq. (^includes contributions coming from two 
distinct processes, P = P^-l-rnr: emission of a photon by 
NP-dipole system with the rate P^, and energy transfer 
to electronic excitations in the NP with the rate Tnr- If 
the system size is sufficiently small, kr <S^ 1, then Tnr 
can be found in a standard way from the longwave limit 
of the Green dyadic?^ 
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is expansion of Coulomb potential in the presence of NP 
over spherical harmonics Yim{Y), and r> and r< are the 
larger and smaller values of r and r' , respectively. The 



where e(w) is the metal dielectric function. As the dipole 
distance to NP surface, d = r ~ R, decreases, P^ and F^^ 
exhibit significantly different behavior. The distance de- 
pendence of Tr is determined largely by the interference 
between direct and NP-induced terms, but it remains fi- 
nite for d/R <^ 1. In contrast, Tnr sharply increases 
for small d due to dominant contribution of high-Z terms 
in Eq. ([5]). Indeed, using Drude form for bulk metal 
dielectric function, the large I approximation for NP po- 
larizability ([7]) is a'/ ~ W7i,/w^, where Wb and jb are, re- 
spectively, bulk plasmon energy and damping rate. For 
d/R ^ 1, the estimate of Tnr then takes the form 
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where I is restricted by Imax ^ kpR 3> 1, kp being 
the electron Fermi momentum. In the local approxi- 
mation when the spatial dispersion of metal dielectric 
function e(w) is neglected, 7b is determined by elec- 
tron bulk scattering length, l^^ due to mainly electron- 
phonon and impurity scattering, 7b ^ vp/lsc, resulting in 
Tnr oc d~^ distance dependence. The RPA calculations 



non-local corrections at small d due to electron- hole pairs Sws (r) is sp-band contribution. The latter is related to 



generatio: 



.43^ 



change this behavior to T^r oc d .In the 



rest of the paper, we present microscopic calculation of 
the decay rates and fluorescence quantum efficiency. 



III. QUANTUM THREE-REGION MODEL FOR 

NOBLE METAL NANOPARTICLE 

POLARIZABILITIES 

For calculation of NP polarizabilities Eq. (U), we 
adopt quantum three-region model developed by us 
previousl y^^'^^ that combines a quantum-mechanical de- 
scription for sp-hand electrons within TDLDA and phe- 
nomenological treatment of d-electrons with dielectric 
function ed{uj) in the region confined hy Rd < R and of 
outside medium with dielectric constant £„. This model 
represents an extension of two- region model (i.e., without 
outside dielectric) that is known to describe reasonably 
well screening effects for symmetric (closed shell) spher- 
ical clusters with electron numbers ranging from several 
dozens to several thousands »^"— In the longwave approx- 
imation, local fields at point r outside NP are determined 
by quasistatic potential 0(r,w) = (j)o{r) + S(j){r,uj), where 



(^,r) 



dV^^i^. 
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is induced potential. The induced density is decom- 
posed as Sn{r) = Sns{r) + Snd{r) + Sn„i{r) , where Sns{r), 
dnd{r)^ and (5n,„(r) are, respectively, contributions from 
sp-electrons, d-electrons, and surrounding medium. The 
sp-electrons contribution is determined from TDLDA 
equation 



47 



dns{r)= / dVP,(r,r') 0(r') + K'K^Ol^^^r') (10) 



where Ps{r,r') is polarization operator for noninteract- 
ing sp-electrons, V^[n{r')] is the (functional) derivative of 
exchange-correlation potential, ^^ n(r) being the ground- 
state electron density; Ps(r, r') and n{r) are calculated in 
a standard way using Kohn-Sham equations.—^ The sys- 
tem is closed by expressing 0(r) in Eq. (fTO)) via 6ns{r). 
This is accomplished by relating (5nd(r) and (5nm(r) in 
Eq. ® back to </.(r) as^^'" 



e2(5nd(r) = V[xd(r)V^(r)], 
e2,5n„(r)=V[xm(r)V0(r)], 



(11) 



where Xdir,uj) = [(ed(w) - l)/47r] 9{Rd - r) is inter- 
band susceptibility with step function 0{x) enforcing 
the boundary conditions and, correspondingly, Xm{r) — 
[(e™ — l)/47r] 6{r — R) is the susceptibility of surrounding 
medium. Using spherical harmonics expansions (/)(r) — 
Y.i^^'Kr)Yira{i) and 5n{v) = ^, 5n«(r)Fi„(f) as weU 
as the continuity of (/>(r) at r = i? and r = Rd, we ob- 
tain decomposition 0^'^ (r) = Wq ' (r) + Sw^ ' (r) -I- 6wi ' (r) , 

where Wq = r^ /e{r) is potential of external field, Swq (r) 
is contribution from d-band and outside medium, and 



the sp-band induced density as^iS 

6wi'\r)= dr'r''^A^^\r,r')Sni'\r'), (12) 



where the kernel A"'{r, r') is given by 
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with e(r) — e^, 1, and e„i for r in regions (0, i?), [Rd, R), 
and (i?, oo). The functions Bq {r,r') and Bs {r,r') 
are, respectively, the l-th harmonics of direct and image 
Coulomb potential: 
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-la'+^\dB'i\Rd,r')\ - pi{rlRd)^[B^\Rd,r') 

^a\l + l)X,r,B'i\Ry)\, (15) 

where /3/(r/i?) = ^R^dRB''^\R,r) is the derivative of 
Coulomb potential at the boundary, 

j3i{x) = Ix-^-^ d{x -1)-{1 + l)2;'6i(l - x), (16) 

and coefficients A^, A^, 77, and a are given by 

Ed - 1 . e™ - 1 
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Tj=l-l{l + l)a2'+iArfA,„, a = Ra/R. (17) 

The combined d-band and outside medium contribution 
has the form 
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We then obtain a closed equation for 6ns as, 
5n^\r) = I dr'r''^P2\r,r') 



wll\r') + Sw^^\r') (19) 



+ / dr"r"^A('\r' ,r")Sni'\r") + V^{r')Sni'\r') 



where the polarization operator Pi {r,r') and the ex- 
change potential^^ V^{r') are calculated from Kohn-Sham 
eigenfunctions." Finally, NP polarizability ^ can be de- 
rived from the induced potential ^ outside NP using the 
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FIG. 1. (Color online) Calculated radiative (a) and nonradia- 
tive (b) decay rates of a molecule near D — 3.4 nm diameter 
Au NP for normal and parallel dipole orientation are com- 
pared to semiclassical (Gersten-Nitzan) decay rates. 



FIG. 2. (Color online) Same as in Fig.[T] but for D = 1.4 nm 
diameter Au NP. 



IV. NUMERICAL RESULTS AND DISCUSSION 



relation (50(') = e2a;/e™r'+^ Using Eqs. ^ and (ITSt . 
we obtain a; — ctQ + a^ , where 



<^ (w) = /i?2'+i [a^'+i Ad(l - A,„) - A,„] /77, (20) 



and 
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Equations (IT^ ~ (P^ determine self-consistently quantum 
polarizabilities ai{uj) of noble metal NP in dielectric 
medium which, in turn, determine the decay rates ([5]) 
and ^. In the next section, we present the results of 
our numerical calculations of Tr and Tnr and of fluores- 
cence quantum efficiency for a molecule near a small NP. 



Calculations of decay rates ([S]) and ([H]) were carried 
out for a point dipole situated near Au NP in a medium 
with dielectric constant e^ = 1-77. We chose dipole's 
emission frequency to coincide with that of Cy5 fluo- 
rophore, 1.83 eV, which lies at low energy side of Au 
NP SP resonance (1.98 eV), so that quantum size and 
non-local effects we are interested in are not being ob- 
scured by interband absorption. To elucidate the role of 
quantum-size effects, the decay rates F^ and Tnr were cal- 
culated for two spherical Au clusters with electron num- 
bers A^ = 92 and N = 1100 corresponding to NP diam- 
eters D « 1.4 nm and D w 3.4 nm, respectively, and to 
highlight the role of non-local and screening effects the 
results of microscopic calculations are compared to those 
obtained using semiclassical polarizabilities, Eq. ([7]) with 
the Landau damping of plasmons 7s = vp/R included 
in both quantum and classical NP polarizabilities. Both 
normal and parallel dipole orientations were considered 
and multipolar polarizabilities ai{uj) with angular mo- 
menta up to / = 30 were included. The intrinsic decay 



rates, F° = 1.8 x 10^ 



and T° = 1.08 X 10^ 



of 



Cy5 dye were included in the full decay rates, and the 
bulk interband dielectric function^ ed(w) was used for 
determination of three-region model parameters i^^*^ 
In Figs. [1] and [21 we show calculated radiative decay 
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FIG. 3. (Color online) Calculated fluorescence quantum ef- 
ficiencies near D = 3.4 nm diameter Au NP for normal 
and parallel dipole orientation are compared to semiclassi- 
cal (Gersten-Nitzan) quantum efficiencies. Inset: blow-up for 
small molecule-NP distances. 



rate F^ and energy transfer rate Tnr for Z? = 3.4 nm and 
D = 1.4 nm Au NPs, respectively, together with the cor- 
responding semiclassical rates based on Gersten-Nitzan 
(GN) model22. defined by NP polarizabilities Eq. (O with 
size-dependent Landau damping correction included. For 
radiative decay rate F^, the non-local effects are not 
significant except in narrow surface region (denoted by 
dashed line) were NP radius is not well defined due to 
spillover effects. This result is expected because only 
I = 1 excitations contribute to F^ through NP dipole 
polarizability ai{uj) which does not account, due to the 
phase space restrictions, for direct excitation of electron- 
hole pairs. A slightly larger quantum F^ for normal 
dipole orientation is due to the difference in d-band 
and sp-band electron density profiles in the NP surface 
layer; the spillover of conduction electrons (as opposed 
to strongly localized d-band electrons) leads to reduction 
interband screening and hence larger local fields in the 
vicinity of NP boundaryi^"— This effect is stronger for 
small NP due to larger volume fraction of underscreened 
regionii^il'^ For parallel dipole orientation, the overall d- 
dependence of F^ is dominated by destructive interfer- 
ence between direct and scattered waves. ^"^ Note that for 
small NPs, the minimum is reached in the spillover re- 
gion, but moves towards outer region with increasing NP 
size. 

In contrast, the calculated quantum energy transfer 
rate F„,. is significantly larger than its semiclassical coun- 
terpart for either dipole orientation. Close to NP surface, 
Vnr for the larger D — iA nm NP is enhanced by fac- 
tor ^ 10 due to direct excitation of Fermi see electron- 
hole pairs by dipole's Coulomb potential^i^ (see Fig.[T]). 
Note that the role of non-local effects in F^^ consider- 
ably exceeds that of SP Landau damping which is in- 
cluded in both calculations. However, Landau damp- 
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FIG. 4. (Color online) Same as in Fig.O but for D = 1.4 nm 
diameter Au NP. 
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FIG. 5. (Color online) Normalized nonradiative decay rates 
for D — 1.4 nm (a) and D — 3.4 nm diameter Au NP are 
shown together with (d/R)'^ fit. 



ing is more important in small D = lA nm NP, as it 
leads to smaller difference between quantum and semi- 
classical Tnr while increasing the overall magnitude of 
Tnr (see Fig. ^ . This results in the overall larger fluo- 
rescence quantum efficiency for larger NPs, as shown in 
Figs. [3] and ID Although for relatively large NP-dipole 
distances (d > R) quantum efficiencies Q of 3.4 nm and 
1.4. nm NPs are comparable, being determined mainly 
by intrinsic molecular decay rates, for small distances 
{d/R ^1) the larger NP Q is significantly higher (see 



insets in Figs. |3] and |4]). At the same time, the difference 
between quantum and semiclassical calculations of Q is 
also larger for D — 3.4 nm NP, indicating a larger role of 
non-local and screening effects. The latter point is illus- 
trated in Fig. ([5]) that compares the distance dependence 
of quantum Tnr with d^^ behavior near a flat surface 
predicted by RPA calculations. '*'^-'*'* In fact, the raise of 
Tnr near NP surface is significantly slower for either NP 
size, but for the smaller D = lA nm NP the deviation 
from d"'* behavior is considerably larger. 

As a final remark, the theory presented here is valid if 
there is no significant overlap between molecular orbitals 
and electron wave functions in the metal. This condition 
no longer holds when molecule- NP distance becomes sev- 
eral A, i.e., in the spillover region (vertical dashed lines 
in Figs. [l]and[2). In this region, fiuorescence quenching 
is dominated by electron tunneling rather than energy 
transfer. We also did not attempt to quantify the effect 
of metal NP on molecules' (or QDs) internal transitions 
and relaxation times which may play important role at 
very close distances i^^i^^ 



V. CONCLUSIONS 

In summary, were performed microscopic calculations 
of radiative and nonradiative decay rates of a fluorophore 
molecule situated near a small An nanoparticle in di- 
electric medium within quantum three-region model that 
incorporates non-local and quantum-size effects and ac- 
counts for different density profiles of d-band and sp-band 
electrons. We found that, close to metal surface, the en- 
ergy transfer rate from molecular dipole to nanoparti- 
cle is significantly (by order of magnitude) higher than 
that predicted by semiclassical electromagnetic models 
due to direct excitation of electron-hole pairs by dipole's 
Coulomb potential. However, the non-local and finite- 
size effects lead to a considerably slower dependence of 
energy transfer rate on molecule-NP distance than d~^ 
behavior predicted by RPA calculations. 

This work was supported in part by the NSF under 
Grants No. DMR-0906945 and No. HRD-0833178, and 
by the EPSCOR program. 
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